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Abstract 

We exhibit a simple and explicit formula for the metric of an arbitrary static 
spherically symmetric perfect fluid spacetime. This class of metrics depends on one 
freely specifiable monotone non-increasing generating function. We also investigate 
various regularity conditions, and the constraints they impose. Because we never 
make any assumptions as to the nature (or even the existence) of an equation of 
state, this technique is useful in situations where the equation of state is for whatever 
reason uncertain or unknown. 

To illustrate the power of the method we exhibit a new form of the "Goldman- 
I" exact solution. This is a three-parameter closed-form exact solution given in 
terms of algebraic combinations of quadratics. It interpolates between (and thereby 
unifies) at least six other reasonably well-known exact solutions. 
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1 Introduction 

The apparently simple problem of the general relativistic static perfect fluid sphere has 
by now generated hundreds of scientific articles. Good summaries of known results, and 
commentaries regarding the extant literature can be found in the book by Kramer et 
al. M, and in the recent review articles by Delgaty and Lake 0, and Finch and Skea ||. 

One of the more common approaches (certainly not the only approach) is to pick 
some barotropic equation of state p = p(p), pick the central pressure, apply the Tolman- 
Oppenheimer-Volkoff equation, and integrate outwards until one reaches the surface of 
the "star" (assumed to be characterized by the innermost zero-pressure surface p — 0). 
Now there are many physical situations in which one simply does not know the equation 
of state, either because of uncertainties in the basic physics (for example, there are still 
some uncertainties regarding the equation of state for nuclear matter in neutron stars), 
or more prosaically because the chemical composition of the "star" may vary throughout 
its bulk so that it is not meaningful to speak of a single equation of state for the entire 
body.Q 

We therefore decided to see what explicit constraints on the spacetime geometry could 
be deduced directly from the perfect fluid condition, without reference to any particular 
equation of state. To start with, note that by using the coordinate freedom inherent 
in general relativity any static spherically symmetric geometry can be put into a form 
where there are only two independent metric components, typically functions of the radial 
coordinate. The most common such forms are given by Schwarzschild coordinates (area 
coordinates, curvature coordinates) 

ds 2 = -|to(f)| dt 2 + g rr (f) df 2 + r 2 {d6 2 + sm 2 6 d0 2 ), (1.1) 

and isotropic coordinates 



ds 2 = -\g u {r)\ dt 2 + g rr (r) dr 2 + r 2 {d8 2 + sin 2 6 A 2 



1.2) 



Now spherical symmetry by itself automatically implies that once one calculates the Ein- 
stein tensor and goes to an orthonormal framed] 

If the geometry is to represent a perfect fluid then in addition we demand pressure isotropy 

Gff = Ggg = Gfa. (1.4) 

This places a single differential constraint on the metric components, and so we expect the 
class of metrics representing a perfect fluid geometry to have only one freely specifiable 

1 A somewhat different approach, explored by Baumgarte and Rendall Q, consists of specifying a non- 
negative but otherwise arbitrary density profile and then integrating the Tolman-Oppenheimer-Volkoff 
equation (without assuming any equation of state) to determine the pressure profile. Such a procedure 
generates geometries more general than the more standard approach sketched above. 

2 Hatted indices will always be used to denote the orthonormal frame attached to a particular coor- 
dinate system. 
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metric component - - more precisely, we expect there to be a single freely specifiable 
generating function, call it z(r), that should characterize the entire class of metrics 



for static perfect fluid spheres. Since the pressure isotropy condition involves derivatives 
of the metric components, we expect the metric g[z(r)] to be some functional of the gener- 
ating function z(r), unavoidably involving derivatives and integrations. These comments 
are of course quite standard and in some form or another implicitly underlie all extant 
static spherically symmetric perfect fluid solutions. The novelty in the current article lies 
in the fact that we will make this implicit procedure explicit and thereby will be able to 
exhibit the most general form of the metric for static spherically symmetric perfect fluid 
spacetimes. That is, we are seeking an explicit closed-form (algebraic-integro-differential) 
solution to the pressure isotropy condition. 

We report that an explicit and relatively simple characterization of this type does 
in fact exist. It involves a single derivative, some algebraic manipulations (of which 
the worst is taking a square root) and an explicit integration. The technique can be 
viewed as a simple algorithm for constructing all static spherically symmetric perfect 
fluid geometries.^ We also discuss the restrictions that must be placed on the generating 
function in order to get "physically reasonable" geometries. 

Finally we present a few specific examples where we demonstrate how various well- 
known solutions fit into our scheme. We exhibit a particularly striking three-parameter 
perfect-fluid solution, given in closed form in terms of algebraic functions. The solution 
is presented in a new manner, and is with hindsight equivalent to the Goldman-I solution 
(Gold-I solution in the Delgaty-Lake classification), which we show is in turn equiva- 
lent to the Glass-Goldman solution (G-G solution). Furthermore, in various regions of 
parameter space the general solution reduces to at least six different previously derived 
solutions. In particular our solution includes three two-parameter sub-solutions: the inte- 
rior Schwarzschild solution, the Stewart solution, and (in the Delgaty-Lake classification) 
the Kuch5 XIII solution. It also contains, as one-parameter branches, the Einstein, de 
Sitter, and anti-de Sitter solutions. We do not claim this list is exhaustive. 

2 Perfect fluid spheres 

Consider a spherically symmetric static spacetime geometry. Without loss of generality 
we know we can put it into isotropic coordinates 



Our first key result can be phrased as a simple theorem: 

3 In particular our technique provides a way of algorithmically generating all possible Baumgarte — 







Rendall configurations. Q 
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Theorem I 

Pick an arbitrary non- increasing function z(r) [that is: a suitably smooth function z(r) 
with z'(r) < 0], introduce a dummy integration variable f, and formally construct the 
metric 



ds 2 



exp { ±2 



+ exp < -2 



V- 



■rz' 



1 — zr 2 



df dt 2 



r ±y/—rz' - 2rz 



zr 2 



dr 



dr 2 + r 2 (d# 2 + sin 2 # 



• (2-2) 



Then this metric is guaranteed to be real [by the non-increasing property of z(r)] and 
always describes a static spherically symmetric distribution of perfect fluid matter. Con- 
versely, the spacetime metric generated by any static spherically symmetric distribution 
of perfect fluid matter can be put into this form for some suitable non-increasing z(r). 

Proof =>- 

By explicit computation 



Gff — Gnp — uj 



[zr 



4y 



g rr (1 — zr 2 ) 2 r 3 



(2-3) 



The computations have been carried out and cross checked using a combination of pencil 
and paper, the CARTAN || package under Mathematical and the standard distribution 
of Maple.Q Invoking the Einstein equations this purely geometric statement ( |2.3| ) implies 



V 



[zr 



4y 



87rG Newton g„ (1 - zr 2 ) 2 r 3 ' 



(2.4) 



QED. 



Proof •<= 

Suppose, on the other hand, we start with a static spherically symmetric perfect fluid. 
Without loss of generality we can put the metric in isotropic coordinates and choose the 
coefficients to be 

ds 2 = - exp{-2<^(r)} dt 2 + exp{+2<^(r) + 4^(r)} [dr 2 + r 2 (dfl 2 + sin 2 d0 2 )j . (2.5) 



Then 



1 

Qrr 



2(p » + Atfj" + ^') 2 + 4(^) 2 + W + ^ + — 



http: / /www, wolfram, com 



http: / /www. maplesoft.com 



(2.6) 



1 

Qrr 



4(^') 2 + ^ - MY 
r 
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(2.7) 
(2.8) 



(^ + < - 2(^') 2 - 







(2.9) 



which is easily solved algebraically (for the derivative if') 

<p' = ±y/2(tp') 2 + (ip')/r - ip". (2.10) 
We could satisfy this equation by picking a "generating function" 0(r) and setting 



9(r) 



V»'(r) 

(p'(r) = ±^/20(r) 2 + 6(r)/r - 9'(r). 



(2.11) 
(2.12) 



But with this particular choice of generating function it is difficult to guarantee the reality 
of the resulting metric. Instead we find it more useful to make the algebraic definition 



ip'(r) 



z(r)r 



1 — z(r)r 2 



that is 



z(r) 



r[l + rip'(r)] 



(2.13) 



With this definition for the generating function z(r) it is now a simple matter to verify 
that the isotropy condition ( 2.10Q is equivalent to 



yf—rz' 



(2.14) 



Integrating and substituting, we get the form of the metric given in the statement of the 
theorem, with now a very simple condition on z(r) [the non- increasing condition] being 
sufficient to guarantee reality of the metric. 
QED. 



Aside: 

We also mention, because it is relatively simple, that for this entire class of metrics 

. V 1 —rz' , , 

P=Tj- a (p + P ■ 2.15 

1 — zr z 

The choice of sign for the square root will be fixed later on, once we demand positivity of 
density at the origin. In contrast, we note that the corresponding formula for G% is quite 
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messy. It is better, but still less than ideal, to consider G% + 3 G r ? = 2i2g which can be 
cast into any of the equivalent forms 



Gff -f- 3 Gf 



5z' + 3r 2 zz' + 2r 3 (z') 2 + r(l - zr 2 )z" 

g rr y/—rz' (1 — zr 2 ) 2 
Az' + {I - zr 2 Y\r(l - zr 2 )- 2 z'}' 

Qrr V~ rz ' (1 — Zr 2 ) 2 

2 



g„ (1 — zr 2 ) r 2 



(r 2 y/—rz')' + y—rz' 



[zr 



4y 



2 

± — 

Qrr 

2 

± — 

Qrr 
Qrr 



(V- 



■rz 



2V- 



■rz' 



1 — zr 2 



r(l — zr 



2\2 



1 — zr 2 
1 



■rz' 



1 — zr' 



V '—rz' 
1 — ^r 2 



r(l — ^r 2 ) 2 
1 I r 2 \/—rz'^ ^ 2zr\/—rz' 



1 — zr 1 



(1 — zr 2 ) 2 



(2.16) 
(2.17) 
(2.18) 
(2.19) 

(2.20) 
(2.21) 



We agree that none of these formulae are stunningly pleasant, but despite considerable 
effort this is (in general) the best we have been able to do. 
Once we apply the Einstein equations 



p + 3p 



87T(jrNewton Q r 



1 / r 2 v '—rz'\' 2zrJ — rz' 
1 1 + 



1 — zr 2 



1 — zr 



2\2 



(2.22) 



Comment: 

We have not used any equation of state anywhere in the derivation. Furthermore, we 
have not yet applied any regularity conditions to the metric — so far it could represent a 
perfect fluid sphere such as a star, or a completely liquid planet (e.g., Jupiter?). It could 
also represent the fluid portion of a mostly liquid planet surrounding a solid core (e.g., 
Saturn), or a black hole surrounded by a spherically symmetric perfect fluid halo (an 
example of a so-called "dirty black hole" or "hairy black hole"), or even a traversable 
wormhole 0, H supported by "exotic" perfect fluid. 

We also wish to contrast this explicit formula for the metric [equation ( |2.2| )1 with the 
more traditional implicit formulation of the problem. (Typically along the lines of "solve 
a certain differential equation for one of the metric components and implicitly substitute 
the result back into the metric ansatz".) 



3 Regularity conditions 

We now investigate the effect of placing various regularity conditions on the geometry 
and the fluid. 
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3.1 Regularity of the geometry at the origin 

If we focus on what is perhaps the astrophysically most interesting case, that of a perfect 
fluid star (or a completely liquid planet), then we want to impose some regularity con- 
ditions at the centre. At a minimum we want the geometry to be regular, which at the 
most elementary level requires 

g tt (r = 0) = finite; g' u (r = 0) = 0; (3.1) 

and 

g rr (r = 0) = finite; g' rr (r = 0) = 0. (3.2) 

We can then without loss of generality rescale r to set g rr ir = 0) = 1 (this only works 
because we are using isotropic coordinates); it is convenient to not rescale g tt (r = 0). 
Then these geometric regularity conditions can be satisfied by: (1) specifying the lower 
limit of integration to be the origin; then (2) setting the integration constants by defining 



ds 2 = - exp{-20(O)} exp J ±2 f /^, dfi dt 2 

\ Jo 1 — zr 2 1 



+ exp { -2 / ± ^ rT7 2fz df \ \dr 2 + r 2 (d6 2 + sin 2 9 d0 2 )l ; (3.3) 
o 1 — zr 2 1 - 



and finally (3) demanding that both 



-> and ^0 as r -> 0. (3.4) 

1 — zr 2 1 — zr 2 

This requires both z(r) and z'(r) to be finite at the origin. 
3.2 Finiteness of central pressure and density 

The central pressure, derived from equation ( |2.4| ) using the condition of geometric regu- 
larity at the origin, is 

Pc = 1 1 , (3.5) 

o7TLt Newton 

which gives no additional constraint beyond regularity of the geometry itself. Indeed 

Z(0) = 27rG'Ncwton Pc- (3.6) 

On the other hand, by considering p + 3p as one nears the origin we can derive addi- 
tional constraints. Suppose we expand z(r) in a power series 

z ( r ) = z (0) + z '(p) r + ^z"{0)r 2 + 0(r 3 ). (3.7) 

Then, evaluating the numerator and denominator of equation ( |2.16|) separately 

|, + 3p](r) =a/ r jm±M&±^l . (3.8) 
87rG N cwton J- rz >(0) - z"(0)r 2 (l + 0(r)) 



Spacetime geometry of static Quid spheres 8 

So if the central value of p + 3p is to be finite we must have 

z'(0) = 0, (3.9) 

in which case 

\P + SplW = * , ^ + „ = + 0(r). (3.10) 

Thus if both pressure and central density are to be finite we must have 

z(r) « ^(0) + ^"(0)r 2 + 0(r 3 ) (3.11) 

with 

±6 



Pe + 3Pc=^ V^'W- ( 3 - 12 ) 

So in summary, finiteness of central pressure and density implies 

^(0)= finite, (3.13) 

while 

*'(0) = 0, (3.14) 

and 

z"(0) = J 8nG ^ 2 [ Pc + 3Pc ] 2 . (3.15) 

3.3 Positivity of central pressure and density 

For the central pressure to be positive we additionally require 

z(0) > 0. (3.16) 

For the central density to be positive we need, first, to take the positive root in the 
expression for p + 3p. This implies a specific choice for the ± throughout the entire body 
of the "star" . That is, we must take 



ds 2 = - exp{-20(O)} exp |2 jT y^^df j dt 2 

+ exp(-2 v/Z ^ 7 ~ 2fZ dfl [dr 2 + r 2 (dfl 2 + sin 2 #d0 2 )l. (3.17) 

[ Jo 1 — zr 2 J L J 

Second, we must demand 

y/-z"(0) > 2z(0). (3.18) 

That is 

z"(0) < -4z(0) 2 . (3.19) 
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3.4 Positivity of pressure (and density?) 

Enforcing positivity of pressure is easy: the pressure is proportional to (zr 4 )' multiplied 
by quantities that are guaranteed positive. So the pressure will be positive as long as 

(zr 4 )' > 0. (3.20) 

The innermost surface at which (zr 4 )' = is defined to be the surface of the star, denoted 

'"surface • 

The positivity of pressure throughout the star then implies 

z{r) + 4—^ > 0; z{r) > -4—^ > 0. (3.21) 

Thus the function z(r) must be positive at least as far out as the surface of the star. 

Guaranteeing positivity of density throughout the star is much more difficult to achieve: 
Mathematically this is because the density (in contrast to the pressure) depends on second 
derivatives z"(r) of the generating function z(r). Physically this arises because we have 
not specified any equation of state. Because of this it should not be too surprising that 
we cannot say everything about the "star" — the surprise perhaps is just how much we 
can do without an equation of state. 



3.5 Monotonic decrease of pressure and density? 

Similarly, guaranteeing a monotone decrease of pressure and density throughout the star 
is generically difficult to achieve. On the other hand, what can be done very easily is to 
derive conditions for the gravitational potential to be monotone increasing as we move 
from the center (that is, to keep the local gravitational force pointing downwards). Given 
the sign choice made to keep the central density positive, we need now only add the 
condition 

z{r) r 2 < 1 (3.22) 

throughout the interior of the star. (This condition also prevents singularities in the 
integration used to define the metric). 

But given our sign choice for the root, we already know 

P' = -T^(P + P)- (3.23) 
1 — zr z 

Thus, assuming monotonicity for the gravitational potential implies 

sign(p') = -sign(p + p). (3.24) 

So under these assumptions the pressure will be monotone decreasing if and only if the 
null energy condition (NEC) is satisfied 0. 
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To calculate the total mass out to radius r we use the Hernandez-Misner mass formula [11 



Mr) = l -R m W = ^-^ 9oe". (3.25) 



This formula is valid for any spherically symmetric spacetime in any spherically symmetric 
coordinate system [HI]. Applied to the metric ( 3.17| ) in isotropic coordinates one obtains 



I — i 2y/-r z'{l + z{r) r 2 ) - 4r z(r) + z'{r) r 2 j 

m(r) = ^(r) r j 2(1 - z(r)r»)2 J " ( } 

The surface of the star is located at r surfaco , where (zr 4 )' = 0. In particular z'(r sur{ace ) = 
— 4^ sur f ace /r sur face- Then a brief computation shows that the total mass of the "star" is 

{2 •\/'2surface ( 1 ^surface '"surface) 4^ sur f ace '"surface 1 
(1 — ^surface r sur face) 2 J 

(3.27) 

Both numerator and denominator can be factorized to yield 



'Ti I, '"surface J V 9O0v surface ) '"surface 



/ \ / / 7 2 '"surface y/ ^surface t-, 2 r S urfacc yj ^surface 

771 V surface,) y 966 v surface ) 7Z , : -^surface 7T - T / = vj - 

^ T y ^surface '"surface ) \ J- T y ^surface ''surface ) 

(3.28) 

Here -Rsurface = \J 9ee ('"surface) is radius of the "star" in curvature coordinates (Schwarzschild 
coordinates). Since each of these factors is manifestly positive, so is the total mass. 
The "compactness" can be defined by 

^ 2 "^('"surface) 4 r SU rface y/ ^surface ^ 2,Q) 

-^surface ~ yj ^surface '"surface ) 

Since 4x/(l + x) 2 < 1, we have x G [0, 1] and the very sensible result 

'"('"surface) < ^surface /2. (3.30) 

That is: 

-^Schwarzschild 2 "^('"surface) 2 (^Newton -^physical ^ -^surface- (3.31) 

These observations are not enough to guarantee that the density is everywhere positive, 
but they do place powerful constraints on its behaviour. 

3.7 Volume-averaged strong energy condition 

Another simple constraint on the stress-energy distribution for a static perfect fluid sphere 
that can be extracted without specifying an equation of state is a certain type of "weighted 
volume average" of the strong energy condition. Specifically, consider 



[ g„ {p + 3p} r 2 dr = — — [ r 2 dr 

J 47T(jr Newton J 



1 / r 2 y/ — rz'\ 2zry/—rz' 
r 2 \ 1 — zr 2 j (1 — zr 2 ) 2 



(3.32) 
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We have already chosen the positive sign for the root, since we want the central density to 
be positive. Note that this is not the "proper volume average" (which would correspond 



to fg, 



3/2 



dr) but is instead weighted by a compensating factor g rr 1//2 , chosen to 



simplify the mathematical analysis. (The occurrence and usefulness of this and related 
weighted volume averages is quite common in spherically symmetric systems.) 

Then, using the regularity conditions we have already deduced for the origin, for any 
value of r* we deduce 



g rr {p + 3p} r 2 dr 



1 



4ttG 



Newton 



r 2 \J—rz' 



1 



zr^ 



2 / r 6 dr 



Zy—rz' 



zr 



2\2 



(3.33) 



The first term is non-negative by the non- increasing property of z(r), which was required 
just to enforce reality of the metric (plus the constraint z(r)r 2 < 1, which was imposed 
to keep the local gravitational force pointing downwards). The second term is positive 
definite for the same reason, so we have, for all r* 



9rr {p + 3p} r 2 dr > 0. 



(3.34) 



This is not the strong energy condition (SEC, p + 3p > 0) itself, but is at least a weighted 
volume average thereof [f|. This implies that the energy density is not permitted to 
become too violently negative. 



3.8 Subluminal speed of sound? 

It is traditional to compute the quantity 





dp 


[dpi 


(8 


fluid " dr 


dr 



(3.35) 



and to demand that this be less than or equal to c 2 , on the grounds that this quantity 
is alleged to represent the physical speed of sound (which certainly should be sublumi- 
nal). This assertion is dangerously misleading, and cannot be justified without significant 
additional technical assumptions above and beyond those that have so far been made. 
Specifically, let us assume that the fluid is described by some equation of state 



p = p(p,x). 



(3.36) 



Here X stands for some collection of variables characterizing the fluid, possibly chemical 
concentrations, entropy density, temperature, or the like. Then 



dp 
dr 



dp 
dp 



x 



dp dp 
oV + dX 



dX_ 

dr 



Thus 






dp 


+ 

X 


dp 


dX 


[dpi 




fluid " 9 P 


dX 


dr 

p 


dr 



1 -1 



(3.37) 



(3.38) 
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That is 



dp" 



c 2 s {X) + 



fluid 



dp 
dX 



'dXT 



(3.39) 



fluid 



In other words (dp/dp)fl ui d can be related to the (constant X) speed of sound c 8 (X) if 
and only if you add extremely powerful additional assumptions. (Such as dp/dX = 0, 
implying an a priori exactly barotropic equation of state. Or dX/dr = 0, implying 
for instance either thorough mixing of the entire fluid mass, an adiabatic star, or an 
isothermal star.) Without such additional assumptions no particular conclusion regarding 
the relationship between (dp/dp)fl U id and the physical speed of sound can be drawn fT0| . 
Given our philosophy in this article (we wish to see what can be deduced without making 
assumptions about the equation of state) such assumptions would be completely opposite 
to our purpose, and so we do not seek to impose the condition dp/ dp < c 2 . 

3.9 Summary 

We can summarize the essential core of these regularity conditions in the following theo- 
rem: 



Theorem II 

Let z(r) be a positive non-increasing function (z'(r) < 0) such that: 

1. z(0) is finite; 

2. z'(0) = 0; 

3. z"{0) < -4z(0) 2 ; 

4. z(r) < 1/r 2 ; 

and consider the metric (guaranteed to be real) 



ds 2 



exp{-20(O)} exp<{2 / / ^ dfldt 2 



o 1 — zr 



f?2 



f n f r V-rz' - 2rz 

+ exp < -2 / — — — dr 

jo 1 — zr 1 



dr 2 + r 2 (d# 2 + sin 2 # 



(3.40) 



Then this metric represents a static perfect fluid sphere with: 

1. regular geometry at the origin; 

2. finite and positive pressure and density at the origin; 



3. a local gravitational field that always points downward. 
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Conversely, any static perfect fluid sphere satisfying these last three conditions can be cast 
into the preceding form with a generating function z(r) satisfying the first four conditions. 

Furthermore, under the conditions enunciated above, if the system is additionally 
isolated (so that the pressure drops to zero at some finite radius), then the total mass is 
guaranteed to be both positive and bounded. 

Proof O 

This theorem is just a codification of the most salient of the preceding results. 
QED. 



4 Examples: 

While the metric given in equation (|3.40|) is guaranteed to be a perfect fluid for a very wide 
class of generating functions z(r), it is only for a much more restricted class of generating 
functions that the relevant integrals can be performed in terms of elementary functions. 
We now present several examples where this can be done. 

4.1 Schwarzschild Exterior Geometry 

The Schwarzschild exterior solution corresponds to 

[mil) 2 

A r ) = "1 ! m = ^Newton M physicah (4.1) 

together with choosing the positive sign for the root. A brief computation leads to p — 0, 
p = 0, and the isotropic form of the Schwarzschild exterior metric. 



, 9 l-m (2r)\ , 9 / m^ 4 

d s 2 = - i£—L at 2 + i + — 

\l + m(2r)J V 2r 



\ 4 

J [dr 2 + r 2 (d6 2 + sin 2 9 d(p 2 )\ . (4.2 



Note that the exterior Schwarzschild does not satisfy the regularity requirements for a 
"normal" fluid sphere — in particular r = is not a "point" but instead corresponds (in 
these isotropic coordinates) to a second asymptotically flat region. Because the geometry 
is not regular at the origin we cannot use equation ( |3.40| ), ( |3.17| ), or even ( |3.3| ). Instead 
we must back-track all the way to ( p.2|) . 

4.2 Einstein universe 



The Einstein universe corresponds to 
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with either sign for the root. A brief computation leads to 

12 4 

P = j P = > (4-4) 

87rG Newton R 2 87rG Newton R? 

and the isotropic form of the Einstein metric: 

_ _ dt 2 , dr 2 + r 2 (dfl 2 + sin 2 £d0 2 ) 

Note that if the density is positive the pressure is negative, and vice versa — the Einstein 
universe does not satisfy the regularity requirements for a "normal" fluid sphere. 

4.3 De Sitter 

The De Sitter geometry corresponds to 

1 3i? 2 + r 2 

* (r) = IvWT^ 2 ' (4 ' 6) 

Choose the positive sign for the root. A brief computation yields 

p _ 12 ^_ 12 

87rG Newton R 2 87rG Newton R? 

and the isotropic form of the De Sitter metric: 

2 f l + r 2 /R 2 \ 2 ^ 2 dr 2 + r 2 (dfl 2 + sin 2 fld0 2 ) 
S " ^l-rVi? 2 ^ + (l-r 2 /i? 2 ) 2 ' 1 ' 8J 

4.4 anti-De Sitter 

The anti-De Sitter geometry corresponds to 

, 1 3R 2 -r 2 

= -R 2 W-^ 2 - ^ 

Choose the negative sign for the root in ( |3.3| ). A brief computation yields 

12 12 

57TLr Ncwton It o7TLr Ncwton JX 



and the isotropic form of the anti-De Sitter metric: 

l - r2,R2 Wt 2 1 dr2+ 

l + r 2 /R 2 J (l + r 2 /R 2 ) 2 



^ = (\zfm\e + ^^^fl. (4.11) 
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4.5 The general quadratic ansatz 

Suppose we consider the general quadratic ansatz 

We were led to this ansatz by considering the form of z(r) for the de Sitter and anti- 
de Sitter cases. By construction, for any choice of Ax, Bx, Cx, and Dx we get a real 
metric satisfying the perfect fluid constraint. (Due to rescaling invariance [multiply both 
numerator and denominator by any fixed constant], only three of these four coefficients 
are actually physically meaningful.) If the central pressure is to be positive then we must 
have z(0) = Ax/ Bx > 0. Then we can without loss of generality take 

*-> = ^Tri£ (413) 

Now in order to justify calling the geometry an "exact solution" we need an explicit 
formula for the metric. Inserting this quadratic ansatz into (|3.40|) the integrals can be 
done in closed form. Expressed in terms of A%, B2, and C2 the resulting metric is a rather 
messy combination of quadratics (in r) raised to various real exponents. 

It is much more convenient to introduce new parameters S, R, and n and write 

2n/(2n 2 -l) 

d, 2 = -l ^'J" 1 dt 2 




( 2 n 2 - 2 n + l)/(2n^l) ^ + ^ + ^ Q #2) 



'l±r 2 /S 2 )(l±r 2 /R 2 



(4.14) 



This is a perfect fluid solution for arbitrary S, R, and n; there are additionally two 
independent sign choices that can be made, one associated with each of the parameters 
S and R. For definiteness of presentation we discuss the case ++; but it is trivial to flip 
the signs as required. 

This solution can (after yet another redefinition of parameters) be seen to be equivalent 
to the Goldman I solution fl^| , called Gold-I in the Delgaty-Lake || classification. The 



pressure and density are rational functions of position. A brief computation yields 

Grr G no G 22 



P 



87rGNewton 87rGNewton 87rGNcwton 

4 qx{r) 



(4.15) 



87rG Ncwton (2n 2 - I) 2 R^ g rr (r) (1 + r 2 / S 2 ) 2 (1 + r 2 / R 2 ) 2 
This verifies that it is a perfect fluid solution. Here ?i(r) is the quartic 

qx(r) = [(2n 2 -1)S 2 R 2 (R 2 -2n 2 S 2 )} 

+ [n 2 (R 2 - S 2 ) 2 - 2S 2 R 2 (2n 2 - l) 2 ]r 2 

-\{2n 2 -l)(2n 2 R 2 - S 2 )]r\ (4.16) 
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Similarly 

G% 4:q 2 (r) 
P = 87rG Ncwton = STrGNcwton (2n 2 - l) 2 R*S 4 g rr {r) (1 + r 2 /^) 2 (1 + r 2 /i? 2 ) 2 ( " } 

where q 2 (r) is the quartic 

q 2 (r) = [3(2n 2 - l)S 2 R 2 ([n - 1}R 2 + n[2n - 1}S 2 )} 

+ [n(n - l)(2n - l)(i? 2 - S 2 ) 2 + 6(2n 2 - lfS^V 

+ [3(2n 2 - l)(n[2ra - l]i? 2 + [n- l]S 2 )]r 4 . (4.18) 

A somewhat simpler quartic is obtained if we consider 

G u + 3Gff An(R 2 - S 2 ) q 3 (r) 



p + 3p = 



87rG Ncwton 87rG Ncwton (2n 2 - l) 2 R*S 4 g rr (r) (1 + r 2 /S 2 ) 2 (1 + r 2 /i? 2 ) 2 ' 

(4.19) 

since then 

q 3 (r) = [3(2n 2 - 1)R 2 S 2 } + [(2n 2 + l)(i? 2 - S 2 )]r 2 - [3(2n 2 - l)]r 4 . (4.20) 

In particular, at the center of the star 

1 R 2 - 2n 2 S 2 

Pc ~ 8vrG Ncwton (2n 2 -l)5 2 i? 2 ' ( ' } 



(4.22) 



while 

1 12[(n- l)R 2 + n(2n- 1)S 2 } 
Pc ~ 87rG Ncwton (2n 2 - 1)S 2 R 2 ' 

and 

- 1 12n(i? 2 - S 2 ) 

Pc + dPc ~ 8vrG Ncwton (2n 2 -l)5 2 i? 2 - ^ 

From this we can use the positivity of central pressure and density to constrain the 
parameters. 

We can find the surface of the star by locating the first zero p(r) [or qi(r)\, and then 
use the Hernandez-Misner formula to deduce the total mass. While a closed-form exact 
algebraic mass formula exists, it is too unwieldy to be reproduced here. 

To summarize the situation so far: The generating function technique developed in 
this note has helped us in several ways. It led us to consider the quadratic ansatz, realise 
it was explicitly integrable, and find a simple form for the metric. We shall now show 
that this quadratic ansatz (the Gold-I solution) is also equivalent to the G-G solution, 
and that it furthermore contains many interesting special cases: interior Schwarzschild, 
Stewart, Kuch5 XIII, de Sitter, anti-de Sitter, and Einstein among them. 
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The G-G [13[] geometry (Glass-Goldman), in the form reported by Delgaty-Lake 



17 



is 



ds 2 



'2r 2 - B 2 D - 1 - C 



2B/C 



2r 2 



B 2 D-1 + C / 
1 



dt 2 



'2r 2 -B 2 D -1-C S 



-(B 2 D+2B-1)/C 



B 2 {2 + D) - (B 2 D + l)r 2 + r 4 \2r 2 - B 2 D - 1 + C 



x 



dr 2 + r 2 (dfl 2 + sin 2 fld0 2 



where 



C 



(B 2 D - l) 2 - 8B 2 



(4.24) 



(4.25) 



Thus at first glance it appears to be a two-parameter solution to the perfect fluid field 
equations. There is a subtlety here in the fact that G-G have implicitly chosen their 
r coordinate to be dimensionless, effectively hiding a dimensional parameter in their 
conventions, that is 

»"this paper = ^0 ?GG (4.26) 

with k,q some arbitrary but fixed distance scale. Then translating the B, D, kq variables 
to our notation 



S 2 
R 2 

n(R, S) 



-4 (l + B 2 D + C)/2; 

C)/2; 



-kI(1 + B 2 D 



I R 2 + k 2 q 
\2S 2 + k 2: 



or 



i s 2 + 4 

\\2R 2 + K 2: 



(4.27) 
(4.28) 

(4.29) 



depending on one's choice for the sign of the square root in the definition of C. That is, 
despite appearances the G-G solution is equivalent to the Gold-I solution and is equivalent 
to our general quadratic ansatz. 

From the point of view of ( |4. 14| ) this is a reflection of the fact that that solution is 
scale- covariant under r — > Ar, S — > \S, R — > \R. 



4.7 Schwarzschild Interior Geometry 

The Schwarzschild Interior geometry is a special case of the quadratic ansatz. It corre- 
sponds to taking both sign-choices positive ++, and setting n — 1. It is now easy to check 
that the metric is 



ds 2 



(1 + r 2 /S 2 ) 2 2 dr 2 + r 2 (d# 2 + sin 2 9 



(l+r 2 /R 2 ) 2 
A brief computation yields 



dt z + 



> 1+r 2/ R 2\2 



P 



8ttG 



Newton 



12 
R 2 "' 



V 



8nG 



Newton 



4 R 2 (R 2 -2S 2 )-(2R 2 - S 2 )r 2 
S 2 R 4 1 + r 2 /S 2 ' 



(4.30) 



(4.31) 
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The central pressure is 

4 R 2 -2S 2 /t n . 

Pc = (4-32) 

OVT Or Newton ti O 

The stellar surface is located at 



R 2 — 2S 2 

Surface = R ]l 2R2 _ $2 ' (4-33) 



The total mass is 



2 

m = — R 
27 



(2R 2 - S 2 ) 3 ' 2 (R 2 - 2S 2 ) 3 ' 2 



(R 2 - S 2 ) 3 

We mention that the generating function is 



(4.34) 



_ B?-2S 2 -r 2 
z[r) ~ R2S 2 + (2R 2 -S 2 y [ j 

Note that S — > is a singular limit of the interior Schwarzschild geometry where the 
central pressure goes to infinity; the central core of the stellar model is on the verge of 
becoming a black hole. On the other hand, as S 2 — ► R 2 /2 from below the stellar surface 
moves inward and the star vanishes. 

If we drive S out of this "regular" range and in particular force S 2 —>■ oo then one 
obtains the Einstein universe (see above). Finally, S 2 — > — Rq and R 2 —>■ +Rq (that is, 
the — h sign choice) corresponds to the anti-De Sitter universe with scale factor R (see 
above); while S 2 — > +R 2 ,, R 2 — ► — Rq (that is H — ) corresponds to the De Sitter universe 
with scale factor Rq. 



4.8 Stewart 



To obtain Stewart's geometry [[L4[] we choose the signs to be and pick n — — 1; it is 

also convenient (but not mandatory) to interchange the roles of R and S. When written 
in this form we can see that it is very closely related to the interior Schwarzschild solution. 
It is now easy to check that the metric is 



(l-r 2 /S 2 ) 2 



dt 2 



[1 - r 2 /R 



2\4 



(1 - r 2 /R 2 ) 2 
A brief computation yields 

1 12R\S 2 



r 2 /S 2 



dr 2 + r 2 (d6 2 + sin 2 9 



„2\4 



P 



and 



V 



[S 2 (2S 2 -3R 2 )-(2R 2 -3S 2 )r 2 ] 

87rG Newton S^ 2 {R 2 -r 2 f 

1 AR\S 2 - r 2 Y[S 2 {2R 2 - S 2 ) + (R 2 - 2S 2 )r 2 } 



8nG 



Newton 



S 12 (R 2 



~2\5 



(4.36) 



(4.37) 



(4.38) 
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The central pressure is 

4 2R 2 -S 2 /t , 

Pc = (4-39) 

oVrCr Newton D 

Which implies i? 2 > S* 2 /2. In contrast, the central density is 

12 2S 2 -ZR 2 
P C = ozp; ^W2~ ' ( 4 ' 4 °) 

57T Or Newton -ft D 

which implies S* 2 > 3R 2 /2. Combined this provides a rather tight constraint 

-S 2 < R 2 < -S 2 . (4.41) 
2 3 

The stellar surface is located at 



2R 2 — S 2 

rsurface = ffJ 2 S 2 - R 2 ' 



The total mass is 



2 Q 

m — — b 
27 



(2S 2 - R 2 f' 2 (2R 2 - S 2 f' 2 



R 4 (R 2 - S 2 ) 
We mention that the generating function is 



(4.43) 



2R 2 -S 2 -r 2 

Z{T) R 2 S 2 + (R 2 -2S 2 )r 2 ' 1 ] 

4.9 Kuchowicz: Kuch5 XIII 

To obtain the Kuch5 XIII geometry |15[ we simply let R — ► oo; it is then convenient (but 
not mandatory) to re-label S as R. It is now easy to check that the metric is 

ds 2 = _(1 + rVi? 2 ) 2 ^ 2 " 2 - 1 ^ 2 + (1 + r 2 /jR 2)(2-2n)/(2n3-l) ^2 + ^2 + ^2 q ^ 

(4.45) 

A brief computation yields 

1 4(n - 1) [3(2n 2 - 1) + (2n 2 - n)r 2 / R 2 } (1 + r 2 / i^^n^-i) / '(2^ -1) 

P ~ 87rG* Newton (2n 2 - l) 2 i? 2 " 

(4.46) 

and 

_ 1 4[(2n 2 - 1) + H 2 r 2 /R 2 ] (1 + 7 .2/ jR 2)-2n(2n-l)/(2n"-l) 

P - 8nG Ncwton (2n 2 - 1)2 R2 • ( 44 ) 

The central pressure is 
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Which implies n 2 > 1/2. In contrast, the central density is 

12 n-1 

Pc ~ 87rG Ncwton i? 2 (2n2-l)' (4 ' 49) 

which further implies n > 1. Because of these constraints there is no stellar surface; 
pressure remains positive for all values of r and the solution is actually cosmological. 
(This ultimately can be traced back to the fact that B2 = 0; which means we are dealing 
with a singular solution of our general three-parameter result.) 
We mention: 

Z{r) = R 2 {2n 2 -\) + 2nV (450) 

Also, it is formally possible to replace R 2 — > —R 2 at the cost of reversing the positivity 
conditions (2n 2 < 1; n < 1). 



5 Discussion and Conclusions 

We have explicitly characterized the spacetime metrics corresponding to the class of all 
static spherically symmetric perfect fluid geometries in a relatively straightforward man- 
ner. This observation is useful whenever there is some uncertainty regarding the actual 
equation of state one wishes to use. The first theorem we presented is applicable to all 
static spherically symmetric perfect fluid geometries without further restriction, while the 
second theorem encodes the most important of the regularity conditions that are relevant 
to an isolated static fluid sphere (such as a star). 

Though the formulae we present do involve an integration, it is particularly noteworthy 
that in our representation the metric is explicit. Furthermore it is easy to keep the 
metric real, and particularly easy to find the surface of the "star". Some (but not all) 
of the standard regularity conditions are easy to enforce, and can be interpreted as extra 
restrictions on the class of "generating functions" z(r). 

Throughout this article we chose to work in isotropic coordinates, because we found 
them to be the most useful. (See Glass and Goldman for an earlier, and rather different, 
use of the ideas of isotropic coordinates and generating functions |H§.) The use of isotropic 
coordinates is not a matter of deep principle and we do not rule out the possibility 
that there may still be other (possibly even simpler) representations in other coordinate 



systems. For instance, the recent work of Fodor [[16] in Schwarzschild coordinates is 
particularly intriguing. 

In closing we reiterate that while a tremendous amount is already known concerning 
static spherically symmetric spacetimes (see in particular |], |], U) the particular approach 
adopted in the present article falls well outside any of the standard schemes.^ 

6 Among the approaches we have encountered in the literature, the closest in spirit to the current 
approach is the Baumgarte-Rendall analysis of H . 
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